Current-activity versus local-current fluctuations in driven flow with exclusion 
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We consider fluctuations of steady-state current activity, and of its dynamic counterpart, the local 
current, for the one-dimensional totally asymmetric simple exclusion process. The cumulants of the 
integrated activity behave similarly to those of the local current, except that they do not capture the 
anomalous scaling behavior in the maximal-current phase and at its boundaries. This indicates that 
the systemwide sampling at equal times, characteristic of the instantaneous activity, overshadows 
the subtler effects which come about from non-equal time correlations, and are responsible for 
anomalous scaling. We show that apparently conflicting results concerning asymmetry (skewness) of 
the corresponding distributions can in fact be reconciled, and that (apart from a few well-understood 
exceptional cases) for both activity and local current one has positive skew deep within the low- 
current phase, and negative skew everywhere else. 

PACS numbers: 05.40.-a,02.50.-r,05.70.Fh 



I. INTRODUCTION 

In this paper we consider fluctuations of the steady- 
state current, and of its close relative, the current activ- 
ity, in the one-dimensional totally asymmetric simple ex- 
clusion process (TASEP). This model is among the sim- 
plest in non-equilibrium physics, while at the same time 
exhibiting many non-trivial properties [J-Q. Some rel- 
evant developments in the study of current fluctuations 
are as follows: exact expressions for the diffusion con- 
stant were found for systems with periodic (PBC) [7j and 
open Q boundary conditions (BC); the full probability 
distribution function (PDF) of current fluctuations was 
similarly considered for both PBC [9] and open [10] BC. 
Very recently, a number of new results have been found 
for current fluctuations in systems with open BC (Il1 - [l4j |. 

In a recent publication [l5[ , exact and numerical results 
were given for steady-state current activity fluctuations 
in the one-dimensional TASEP, for both periodic and 
open boundary conditions. By making use of the known 
steady-state operator algebra [3fl , exact expressions were 
derived for the three lowest moments of the activity PDF, 
which fully display their finite-size dependence. All these 
were confirmed to excellent degree of accuracy by numer- 
ical simulations. The results of Ref. [l5| extend and com- 
plement earlier analytic work on the joint distribution of 
current activity and density for the TASEP We 
recall that the current activity (henceforth denominated 
simply activity) is not identical to the standard current, 
although the first moments of the respective distribu- 
tions coincide. As explained below, the former quantity 
is static, in this sense akin to the instantaneous (local or 
global) particle density, while the latter is a dynamic one. 

Many exact results available for current fluctuations 
pertain to the infinite-system limit 0, [l(| EH , although 
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the diffusion constant has been calculated for finite sys- 
tems 0, Finite-size effects have been considered also, 
e.g., in Refs. (TMl- 

Our main purpose here is to exploit the possible 
connections between activity- and current fluctuations. 
Given that the former quantity has proved amenable to 
such detailed description, it is desirable to check whether 
its properties can help explain any relevant aspects of the 
latter. 

In the time evolution of the 1 + 1 dimensional TASEP, 
the particle number n/. at lattice site I can be or 1, and 
the forward hopping of particles is only to an empty ad- 
jacent site. The stochastic character comes from random 
selection of site occupation update [H, [l9j : if site I is 
chosen for update, the instantaneous current across the 
bond from t to I + 1 is given by Jgj+i — ng(l ~ n? + \) . 

With open boundary conditions, the case considered 
here, the additional externally-imposed parameters are: 
the injection (attempt) rate a at the left end, and the 
ejection rate /3 at the right one. The phase diagram in 
a—f} parameter space is known exactly, as well as many 
other steady state properties lis LLZL 120, ■ 

The total (instantaneous) activity A within the system 
is defined as the number of bonds that can facilitate a 
transition of a particle in the immediate future. Thus 
it equals the number of pairs of neighboring sites that 
have a particle to the left and a hole to the right [lfl [l7j • 
For systems with open BC, one usually includes also the 
injection and ejection bonds at the system's ends, though 
these have to be weighted by the respective injection and 
ejection rates, a and /3. 

For an L-site system with open BC (L sites and L + 1 
bonds, including the injection and ejection ones), one has: 

L-l 

A = a (1 - m) + rij (1 - n e+1 ) + f3 n L (Open BC) . 

(1) 

The activity is a snapshot of the system at a given 
moment in its evolution; in this sense, it is as much of a 
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static quantity as, for example, the instantaneous global 
density. By contrast, the current is a dynamic object, as 
it reflects the stochastically-determined particle displace- 
ments which actually take place during a unit time inter- 
val. The investigation of current fluctuations is usually 
carried out by examining the total charge (i.e., number 
of particles) crossing a given bond, during a long time 
interval in the steady state regime [8l-ll4|. 

The (properly normalized) first moments of activity 
and current PDFs coincide. For A as defined in Eq. ([T]) 
one has: 

J = j^j (A) (Open BC) , (2) 

where J is the average steady-state current through any 
bond in the system, and brackets (...) denote ensem- 
ble averages. The above equality can be understood by 
recalling that successive steady-state snapshots (activity 
configurations) are generated via the intervening particle 
hoppings, which constitute realizations of the system's 
current. In our simulations of Ref. [TH we verified that 
this property holds, to within numerical accuracy, in all 
cases investigated there. However, the connection at this 
level is not sufficient to warrant equality of higher mo- 
ments of the PDFs. 

Section [TT] below recalls selected existing results, and 
gives a theoretical background to the concepts used in 
this work. In Section imi our numerical simulations are 
described, and their corresponding results are exhibited. 
In Sec. IIVI we provide a global analysis of the numerical 
results; finally, concluding remarks are made. 

II. THEORY 

In studies of current fluctuations for the TASEP, one 
considers the steady-state current through a specified 
bond connecting sites I and I + 1, henceforth denoted 
by Jg, and the associated integrated charge Qe(t) = 
f*J e (t')dt'. With open BC, the leftmost (injection) bond 
£ = is usually singled out for examination, although the 
results for the PDF of steady-state current fluctuations 
would be the same for any choice of i Q. Thus, here we 
frequently omit bond labels wherever this does not give 
rise to ambiguity. Also, it is usual to remove the linear 
term from the integrated charge, and to consider instead: 

Q(t) = Q(t) - Jt , (3) 

so (Q(t)) = 0. Closed-form expressions are available for 
J as a function of a, /3, and number of sites L [3|; see, 
e.g., Eqs. (29) and (32)-(34) of Ref. OJ. 

For t — > oo, exact expressions for the lowest-order cu- 
mulants, ((Q n )), of the integrated current have been pro- 
posed: for n — 2 and any L, everywhere on the phase 
diagram Q ; for n = 1 — 3 and L — > oo , anywhere except 
for the maximal-current phase a, (3 > 1/2 [l4j|; and (via 
a parametric representation) for any a, /3, and L, essen- 
tially for all n [l3| ■ In the latter, explicit expressions are 



forthcoming only for L — > oo away from the maximal- 
current phase, and for n < 3, any L, at (a, f3) = (1, 1). 
The above results agree with each other, wherever com- 
parison is possible. Furthermore, on the basis of gen- 
eral properties of the associated generating function (see 
e.g. Ref. 01), the assumption is made that all cumulants 
scale linearly with time, so the quantities of interest are 
the E n = ((Q n ))/t, which are frequently referred to as 
cumulants of the local current (l4| . 

A direct elementary illustration of such linear de- 
pendence can be given by considering the variance 
((Q 2 (t))) . With 8Jg(t) = Jg(t) — J, assuming an ex- 
ponential decay of the current-current correlation func- 
tion, (5JfXt')5J e {t")) ~ e-l*'- t "l/ r , and using steady- 
state properties, one has: 

((Q 2 (t))) = [ dt' [ dt"(6Mt')6J e (t")) = 
Jo Jo 

(for i > r) = i / ds(SJ e (0)SJ e (s)) = tr . (4) 
Jo 

An interesting exception to the above has been pointed 
out in Ref. I 111 , where theoretical and numerical argu- 
ments are presented to show that ((Q n )) ~ t ™/ 3 for n > 1 
at the second-order phase boundary between low- and 
maximal-current phases. We shall return to this later. 

As recalled in Eq. (0]), the cumulants ((Q n )), n > 1, 
involve unequal-time correlations @, because the inte- 
grated charge accumulates local current fluctuations over 
time. In contrast, see Eq. (JXJ) , the nontrivial features of 
(instantaneeous) activity statistics arise because it adds 
contributions from all sites in the system, thus it de- 
pends on non-local, equal-time, correlations (l5| . In fact, 
it is because A(t) is the sum of O(L) equal-time stochas- 
tic (albeit not independent) variables that the variance 
and skewness (and, most likely, higher-order moments) 
of its PDF always approach zero for large L [l5j . The 
spatially local character of the ((Q n )), in turn, implies 
less severe L-dependent effects on these quantities; thus 
the E n generally converge towards non-zero values for 

oo [IIiIClI]. 

We wish to consider the relative importance of local 
and non-local (equal and non-equal time) correlations. 
We then define a hybrid quantity, the position-averaged 
(instantaneous) current J(t): 

e=o 

The ensemble average of Jit) coincides with that of the 
normalized activity, see Eqs. (HJ and ([2]), and of course 
with the steady-state current J. The cumulants of its 
integral over time, Q(t) = L J(t') dt' incorporate both 
equal- and unequal-time (as well as local and non-local) 
correlations. For consistency, see Eq. ©, we also sub- 
tract the linear contribution, Jt, from Q(t). 

In order to draw relevant distinctions and similarities 
with Q(t), we also consider the integrated (normalized) 
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activity, A(t) = (L + l)" 1 f*A(t/)dt' (again, with the 
linear contribution Jt subtracted). 

In what follows, we show results of numerical simula- 
tions of J(i) and associated quantities such as the cu- 
mulants of Q(t), as well as those of A(t). in comparison 
with those corresponding to the local current. 



III. NUMERICAL RESULTS 
A. Introduction 

To make contact with previous work, we usually con- 
sidered lattices either with L = 600 sites, as done in 
local-current simulations or with L — 256, for which 
many results on activity statistics are available . It is 
important to take L > 1 not only in order to minimize 
finite-size effects, but also because some features such as 
the non-linear scaling of current cumulants with time, re- 
ferred to above, occur only during time "windows" whose 
width increases as L grows large [Tlj |. 

In our simulations, a time step is defined as a set of 
L sequential update attempts, each of these according 
to the following rules: (1) select a site at random; (2a) 
if the chosen site is the rightmost one and is occupied, 
then (3a) eject the particle from it with probability /3; 
alternatively, (2b) if the site is the leftmost one and is 
empty, then (3b) inject a particle onto it with probability 
a: finally, if neither (2a) nor (2b) is true, (2c) if the site 
is occupied and its neighbor to the right is empty, then 
(3c) move the particle. 

Thus, in the course of one time step, some sites may 
be selected more than once for examination, and some 
may not be examined at all. This corresponds to the 
random-sequential update procedure of Ref. [IH 

For the various sets of a, f3 considered here [because of 
particle-hole duality, we take only a < (3 ] , we found that 
for L = 600 and starting from an initial random config- 
uration of occupied sites, one needs n m = 10, 000 time 
steps for steady-state flow to be fully established, so that 
the ((Q n )) are free from startup effects (for a = j3 = 1 
this remark needs to be qualified, see Section MI CI) . 
Hence we discarded the first nj n configurations when eval- 
uating all quantities discussed in the following. As seen 
below in Figures [TJ [5J and the characteristic decay 
times for the ((Q™)), ((A n )) are about twice that, so 
the respective curves still show some crossover towards 



pure power-law behavior for the first 
depicted. 



10 time steps 



B. a + P = 1 

In Figure Q] we present the evolution of ((-4 2 )) and 
((Q 2 )) against time, for (a,/3) = (0.3,0.7) and (0.5,0.5). 
In both cases the variance initially increases faster than 
linearly with time, but eventually settles to a linear 
increase. Also, both quantities behave identically for 
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Figure 1. (Color online) Double-logarithmic plot of variances 
((.A 2 )), for the integrated activity A(t), and ((Q 2 )} for Q(t), 
both defined below Eq. ([5]), against time. The dashed line is 
oc t. System size L = 600. 



a = 0.3, (3 = 0.7, while for the special case a = /3 = 0.5 
they follow very close trends, the ratio ((A 2 )) / ((Q 2 )) 
starting at around 0.85 and approaching unity as t grows 
(near the upper limit in the Figure, it has reached 0.95). 

One sees that the dominant features of the second- 
order cumulants of Q essentially coincide with those of 
A. As shown below, the quantitative description of n = 3 
cumulants, and of their associate skewness properties, al- 
ways points to their absolute value being rather small; at 
the level of accuracy pursued here the main issue there- 
fore concerns their sign, which generally proves to be a 
robust feature (i.e., independent of whether Q or A is 
the quantity under consideration). Thus, in the follow- 
ing, we usually display only results for the ((<2™)), for 
comparison with the conventional ((Q n )) ■ 

In Figure[5]we show the cumulants ((Q 2 )) and ((Q 2 )) 
against t, for a = 0.3, /3 — 0.7, and for a — (3 — 0.5 . 

Regarding data for ((Q 2 )) in Fig. for a = 0.3, /3 = 
0.7 the expected linear behavior against t is present for 
the whole time interval shown. On the other hand, for 
a = f3 = 0.5 one sees the anomalous scaling with i 2 / 3 , 
referred to above, over a wide time window up to t sa 
3 x 10 4 . Both types of power-law behavior have been 
reported in Ref. [Ly, for the same two sets of values of 
(a,j3). For the a = /3 = 0.5 data, at longer times there 
is a crossover back towards linear scaling. 

Anomalous scaling at a = ft = 0.5 has been explained 
in detail in Ref. It can be understood by recalling 
that this point is on the second-order phase boundary 
separating low- and high-current phases (respectively for 



4 




Figure 2. (Color online) Double-logarithmic plot of variances 
{{Q 2 )), for Q(t) defined in Eq. ©, and ((Q 2 )} for Q(t) defined 
below Eq. (JSJ , against time. The full line is oc t 2//3 ; the dashed 
line is oc t. System size L = 600. 



Figure 3. (Color online) Scaling plot of PDF for the variable 
Q(t) defined below Eq. ©, at {a, /3) = (0.5,0.5). iV sam = 10 6 
steady-state samples were collected at the indicated times 



a < 1/2 or (3 < 1/2, and a, (3 > 1/2). The associated 
diverging correlation length £ brings about a diverging 
relaxation time r ~ £ z (where z = 3/2 is the Kardar- 
Parisi-Zhang exponent known to describe TASEP dy- 
namics 0, 0, E|) which, in turn, governs the current- 
current correlation fluctuations and related quantities. 
For systems of finite size L, this regime holds only for 
t < L z , i.e., the "windows" referred to earlier. In con- 
nection with Eq. Q, it can be seen that in this case, the 
assumption of t 3> r fails over the window of anomalous 
scaling, but is then restored at longer times; hence, the 
observed crossover towards linear behavior. 

Going now to data for ((Q 2 )), one sees that they 
eventually merge with the respective ((Q 2 )) curves. For 
a = j3 = 0.5 such merging happens, of course, at later 
times than the window of anomalous scaling. It was 
shown in Ref. [l5| that the activity PDF is a pure Gaus- 
sian at this point; here, as already remarked in connec- 
tion with Figure [TJ we see that ((<2 2 )) does not capture 
the anomalous scaling behavior there, either. 

Further information can be extracted from the PDFs, 
or full counting statistics, of the integrated currents. 
Data for Q(t) at (a, 13) = 0.5 are shown in Figure [3] as 
a scaling plot. Although the effective scaling power is 
1/2, which is in line with the linear dependence of ((<2 2 )) 
against t seen in Fig. [5J the scaled PDF shows a slight 
negative skew. For comparison, Figure 0] shows the cor- 
responding scaling plot for Q(t) within the anomalous 
scaling window, where a negative skew is more clearly 
present; however, the scaling power is now 1/3, as al- 
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Figure 4. (Color online) Scaling plot of PDF for the variable 
Q(t) defined in Eq. ©, at (a, /3) = (0.5,0.5). iV sam = 10 6 
steady-state samples were collected at the indicated times 



ready remarked in Ref. [Til 

For a proper analysis of the scaled data shown in Fig- 
ures[3]and[31 and comparison with theoretical predictions 
where available, one must keep in mind that the former 
necessarily pertain to finite t and L, while the latter ei- 
ther assume taking the L,t —> oo limit [H, or Q L 
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is kept fixed and t oo. 

Denoting by k the scaling powers referred to, it 
turns out from the scaled variables u = Q/t k \Q/t k ] 
and P(u) = t k P(Q) [t k P{Q)\ that the cumulants of 
the global [ local ] current are given by £„ [ E„ } — 

t nk-l ^ u n^ _ 

So, for the global current (k — 1/2), one can have both 
((it 2 )) and £2 finite for t — > 00, while ((ii 3 )) must vanish 
in the same limit, since £3 must not diverge. By sepa- 
rately analyzing PDF data used in Figure [3] for various 
times 4000 < t < 24000 we found that, as t increases, 
£ 2 = ((u 2 )) increases from 0.0077 to 0.0081, while ((w 3 )) 
goes from -2.7 x 10~ 4 to -1.4 x 10~ 4 . The "diffusion 
constant" £2 can be compared with E2 « (167T-L) -1 / 2 = 
0.0088 ... for i = 256 d|. Although there appears to be 
no a priori reason to assume equality of these quantities, 
the closeness of their values is remarkable. As regards 
((it 3 )), its evolution can be approximately matched by 
a t -1 / 2 dependence, which would be compatible with a 
nonzero limiting value for £3. 

For the local current (k = 1/3), if ((u 2 )) ^ then 
E2 = ((Q 2 )) Jt ~ i -1 / 3 which is in line with existing re- 
sults IE Mill, recalling that at (a, ft) = (0.5,0.5): (i) 
E2(L) oc L~ x l 2 [H, and (ii) the appropriate finite-size 
scaling combination is t/L 3 ^ 2 . Also, E3 should approach 
a finite value if ((it 3 )) is finite. From the full set of data in 
Figure [31 we find ((it 3 )) = —0.034(1) . Separate analysis 
of subsets of data for different times shows small fluctu- 
ations around the average just quoted, with no apparent 
up- or downward trend upon increasing t. Comparison of 
this estimate to the analytic prediction for the thermo- 
dynamic limit, E3 — —1/8 [IE HI shows that, though 
correct in sign, it is still only about 1/4 of the expected 
value. 

At (a, (3) = (0.3,0.7), as evinced in Figured the 
scaling power is k = 1/2 for both Q and Q. Fig- 
ure [5] shows that the scaled PDFs for both quantities are 
nearly indistinguishable, even when one plots L = 256 
data for P(Q) together with L = 600 data for P(Q). 
From fits of pure Gaussians to the scaled PDFs, we get 
£2 ~ E2 — 0.0841(5), in very good agreement with 
the thermodynamic-limit prediction [E [IE HI E2 — 
a(l — a)(l — 2a) = 0.084. Allowing for non-zero skew- 
ness gives E 2 ~ 0.082, with E 3 < 0, and \E 3 \ < 1 x 10~ 3 . 
However, the estimates for the latter quantity exhibit un- 
certainties of order 60% or thereabouts, so they should be 
considered with caution. The theoretical prediction for 
L -> 00 HlHl is E 3 = a(l-a)(l-6a+6a 2 ) = -0.0546 . 
As in the case of (a, ft) = (0.5, 0.5), we get the right sign 
for E3 but a smaller absolute value than predicted. 



C. a = P = l 

Within the maximal current phase a, ft > 1/2, the 
behavior is expected to be similar to that observed at 
(a, ft) = (0.5,0.5), and reported in Section UlIBI At a — 
ft = 1 several simplifications occur 

[El, allowing 
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Figure 5. (Color online) Scaling plot of PDFs for: (i) the 
variable X — Q(t) defined in Eq. © [ system size L — 600], 
and (ii) X = Q(t), defined below Eq. © [L = 256], at 
(a, ft) = (0.3,0.7). The full (red) line is a Gaussian fit to 
data. iV sam = 5 x 10 J steady-state samples were collected at 
the indicated times . 



for simple expressions to be derived for the ((Q n )), so we 
ran simulations at that point. 

Figure [S] shows the time evolution of the variance of 
the several integrated quantities under investigation here. 
Comparing with Figures [TJ and [2J which refer to the same 
system size and time interval, the main difference to the 
behavior exhibited at (a, ft) = (0.5, 0.5) is in the anoma- 
lous scaling of ((Q 2 ))- Our best fit to single-power be- 
havior, encompassing the same time "window" as used 
in that case, i.e., 200 < t < 3 x 10 4 , gives an exponent 
equal to 0.61(1), close to but slightly below the expected 
value, 2/3. The expected crossover to linear behavior, 
and merging with the ((<2 2 )) curve, takes place in the 
same (narrow) time interval as at (a, ft) = (0.5,0.5). 

The apparent discrepancy in the exponent value can be 
solved by examining the scaling plot of the PDF for Q, 
shown in Figure where the scaling power is set to 1/3. 
It can be seen that data collapse is rather good, except 
for the t = 4000 data which visibly stray off, especially 
at the low end of the curve. Indeed, fitting only data 
for 8 x 10 3 < t < 3 x 10 4 in Figure [5] gives an exponent 
0.66(2). So it is the extent of the i 2 / 3 scaling window 
which is shorter here than for (a, ft) = (0.5, 0.5) . Further 
analysis of the scaling plot depicted in Figure [7] follows 
the same lines as that of Figure HJ (i) again, E2 (t) ~ 
t -1 / 3 , with the same interpretation as above; and (ii) the 
value of ((it 3 )) = —0.012(1) compares reasonably well to 
the analytic prediction, lim£_ i . 00 E3 = —0.009 . . . [13]. 

Figure [5] shows the scaling plot for the PDF of <2, using 
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Figure 6. (Color online) At a = /3 = 1: double-logarithmic 
plot of variances {(Q 2 )), for Q(t) defined in Eq. ©, ((Q 2 )} for 
Q(t), and ((.4 2 )) for the integrated activity A(t), the latter 
two defined below Eq. ([5}, against time. The full line is oc 
t ' 61 : the dashed line is oc t. 



Figure 7. (Color online) Scaling plot of PDF for the variable 
Q(t) defined in Eq. ©, at (a,/3) = (1,1). iV sam = 5 x 10 5 
steady-state samples were collected at the indicated times. 



1/2 as the scaling power. The quality of data collapse is 
not as good as for points at the border, or outside, of 
the maximal-current phase (see, respectively, Figures [3] 
and[S]). However, it is possible to find a rather stable esti- 
mate for the "diffusion constant", £2 = 0.0045(1). Again, 
this compares favorably with E2 — (3v2tt/64) L~ 1//2 = 
0.004797... for L = 600 [I3|. The third cumulant is 
negative but very close to zero, £3 « — 1 x 10~ 4 . 



IV. DISCUSSION AND CONCLUSIONS 

We have set out to compare the features exhibited by 
fluctuations of activity with those of (local) current for 
the TASEP with open BCs. 

We define C„ = ((A n )), with the systemwide (instanta- 
neous) activity A defined in Eq. [TJ for compatibility with 
the notation of Ref. [T^. So C„ = C n /(L + 1)™ is its nor- 
malized counterpart. As shown in Ref. H5L C2 vanishes 
with system size as L^ 1 for any (a,/3). The standard 

skewness S, defined as S = C^/C^ 2 = C^/C^ 2 (H|, van- 
ishes as L~ x (so C3 ~ L - ( x+3 / 2 )), usually with x = 1/2 
(x = 3/2) in the low (high)-current phase, except for 
some special lines or points on the phase diagram (see 
Figure 8 of Ref. El). As remarked in Section HH C2 and 
C3 vanish with increasing L because they reflect fluctu- 
ations of the equal-time, nonlocal correlations of 0(V) 
stochastic variables. 

On the other hand, cumulants of the integrated activ- 
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Figure 8. (Color online) Scaling plot of PDF for the variable 
Q(t) defined below Eq. ©, at (a, /3) = (1,1). 7V sam = 10 6 
steady-state samples were collected at the indicated times 



ity A, as well as those of the position-averaged current 
J , behave similarly (albeit with different decay times on 
the approach to steady state, see Figures [5] and [5]) to 
those of the local current, except that they do not cap- 
ture the anomalous scaling which occurs for a, (3 > 1/2. 
This indicates that the systemwide sampling at equal 



7 



times overshadows the subtler effects which come about 
from non-equal time correlations, within the high-current 
phase and along its border. Anomalous scaling can be un- 
veiled only when the local current is considered, in which 
case systemwide, equal-time sampling is absent, and the 
system-size dependence is much less relevant overall (al- 
though it still shows up in the width of the scaling "win- 
dow", see Ref. I 111 and Figures [5] and [5] above) . 

Going back to the skewness of distributions, the sign 

of S is positive for min (a, ft) < ~ - *J = 0.27639 
and negative otherwise, with the following exceptions: 
(a, 13) = (0.5,0.5) where S = 0, and (a, (3) = (1,1), 
(1,0.5), and (0.5,1) where S > [Hj]. At min(a,/3) = 
0.27639 . . . , S > but depends on L~ 3 / 2 , as opposed to 
the L" 1 / 2 form which holds generally in the low-current 
phase. All such exceptions are accidental, i.e., they occur 
because of symmetries or cancellations which are valid 
only for those specific values of (a, (3) [l5|. Further- 
more, it has been shown that S > is to be expected 
on physical grounds for a <C 0.5 in the low-current, low- 
density phase (see Figure 3 of Ref. HH and accompanying 
remarks) . 

Concerning local-current statistics, for the low-current 
phase Refs. [lj and [Til predict that, with 7 = min (a, /3), 
E3 = 7(1 — 7)(1 — 67 + 67 2 ) which is thus zero for 70 = 

i — ^5 =0.2113... , and positive (negative) for 7 < 70 
(7 > 70). Ref. [H evaluates E 3 < at (a,/3) = (1,1), 



as recalled in Section UlI C[ and gives arguments showing 
that the dominant behavior of the cumulants should be 
the same as that at (a, j3) — (1, 1) everywhere within the 
maximal-current phase. 

In conclusion, once accidental exceptions are under- 
stood as such, both S for activity and E3 for the local 
current tell essentially the same story: fluctuation distri- 
butions have positive skew for low injection (or ejection) 
rates, and negative skew otherwise. The borderline be- 
tween the two types of behavior lies deep within the low- 
current phase, and its relationship to the actual (second- 
order) transition to the maximal-current phase it at best 
that of a precursor. 

As a final remark, we emphasize that accurate numeri- 
cal checks of the theoretical predictions for E3 would need 
much longer simulations than those reported here, since 
evaluation of this quantity strongly depends on a proper 
description of the tails of the local-current PDF. How- 
ever, at the level of accuracy pursued here, the present 
results fulfil our goal of providing a comparative anal- 
ysis of the main features of activity- and local-current 
fluctuations. 
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